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ABSTRACT 

Context. We examine the evolution of low- mass star and brown dwarf eclipsing binaries. These objects are rapid rotators 
and are believed to shelter large magnetic fields. 

Aims. We suggest that reduced convective efficiency, due to fast rotation and large field strengths, and/or to magnetic 
spot coverage of the radiating surface significantly affect their evolution, leading to a reduced heat flux and thus larger 
radii and cooler effective temperatures than for regular objects. 

Methods. We have considered such processes in our evolutionary calculations, using a phenomenological approach. 
Results. This yields mass-radius and effective temperature-radius relationships in agreement with the observations. We 
also reproduce the effective temperature ratio and the radii of the two components of the recently discovered puzzling 
eclipsing brown dwarf system. 

Conclusions. These calculations show that fast rotation and/or magnetic activity may significantly affect the evolution 
of eclipsing binaries and that the mechanical and thermal properties of these objects depart from the ones of non-active 
low-mass objects. We find that, for internal field strengths compatible with the observed surface value of a few kiloGauss, 
convection can be severely inhibited. The onset of a central radiative zone for rapidly rotating active low-mass stars 
might thus occur below the usual ~ 0.35 Mq limit. 
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1. Introduction 



Low-mass stars (LMS), i.e., M-type stars in the present 
context, represent an overwhelming fraction of the Galactic 
stellar population (jChabrie r 2003) . Observational determi- 
nation of their mass-radius relationship provides a strin- 
gent testing of the theoretical description of their structure 
and evolution. Although observations of this relationship 
agree well with the theoretical predictions in a large num- 
ber of cases ( [Segransan et al. 2003 IChabrier et al. 20"05l) . 
the radius determinations of eclipsing binaries (EBs) 
depart significantly ( ^10-15%) from thes e predictions 
([Torres fc Ribas 2002J IChabrier et al. 2005^ . These dis- 
crepancies are unlikely to be due to inaccurate equation 
of state or opacities, whose treatment is relatively well 
mastered for to^O.4M0 and Toff > 4000 K (Chabrier & 
Baraffe 2000, CBOO). Furthermore, active LMS are ob- 
served to be redder, i.e. cooler, than the other LMS 
( [Hawley et al. 19961 ). These particular behaviours pose a 
challenge to theorists. On the other hand, the majority of 
LMS, notably EBs, are fast rotators and exhibit strong per- 
sistent coronal (X-ray) and chromospheric (H^) activity or 
large flares (jCizis et al. 2002 1. indicating the presence of 
a strong magnetic field (jPonati et al. 2006| . In this Letter, 
we suggest that the discrepancy between the observed mass- 
radius relationship of EBs and of other LMS does not stem 
from inaccurate equation-of-state or opacity problems, but 
is due to a rotation or magnetic field induced reduction of 
the efficiency of large-scale thermal convection in their inte- 
rior, leading to less efficient heat transport. We show that 
the reduction of the star's radiating surface due to mag- 
netic spot coverage also yields a smaller Toff and a larger 
radius. 



2. Effect of rotation and magnetic field 

The EBs and, more generally, active LMS are fast rotators 
with rotation periods of P < 3 days and rotation velocities 

of Vrot >10kms-i (|Delfosse et al. 19981 iReid et al. 20021 
|Mohanty fc Basri 2003[ ). Typical convective time scales in 



their interior are tr 



Hp 



10" 



40^-10^ s, where 



(lO^-lO'i ) 

Hpr^Ri, is the pressure scale height ()CBOO|) . This yields 
^<10-2. 



Rossby numbers Ro= 



At a very small Rossby 



number, the Proudman- Taylor theorem enforces fluid mo- 
tions to columnar or sheet-like structures with a charac- 
teristic length scale perpendicular to the rotation axis that 
is much smaller than the one along the rotation axis (of 
the order of i?*). Although the conditions of this theorem 
do not exactly apply for spherical objects with finite vis- 
cosity, rotation-dominated convective motions are severely 
affected by rotation, leading to highly anisotropic patterns 
with large-scale fluid motions being confined to columns 
along the rotation axis and motions in the other directions 
being strongly inhibited ( Zhang fc Jones 1997[ ). As a whole, 
this reduces the characteristic length scale along these di- 
rections and the mean convective velocity amplitudes, thus 
the efficiency of thermal convection to transport the inter- 
nal heat flux. 

On the other hand, fully ionised interiors of LMS are 
excellent electrical conductors. The interaction between 
the magnetic field and the fluid motions induces an elec- 
trical current and a Lorentz force, which in turn affects 
the motions. In the presence of a magnetic field, the 
anisotropy of rotation-dominated convection is reduced and 
the large-scale flows reach the so-called MAC (Magnetic, 
Archimedean, Coriolis) balance between the buoyancy. 
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Coriolis, and Lorentz forces (jStarchenko fc Jones 2002| . 
Equipartition between the buoyancy force, ~ pgS{AT/T) 
(where ^=~{§^^)p and AT is the temperature excess over 
a pressure scale height), and the Coriohs force, ~ 2pQvconv 
(where $7 — is the star angular velocity), yields 



(1) 



for the typical average large-scale velocity. For f2 ~ 10"'' 
s""'^, typical values of the various quantities over the LMS 
domain give Vconv ~10-100 cms~^, about a factor of 10 
lower than the predictions of the usual mixing length the- 
ory. 

There is presently no clear understanding of the inter- 
action between convection and magnetic field under stellar 
interior conditions. We can, however, estimate the condi- 
tions for generation of a dynamo in LMS or BD interi- 
ors. The microscopic magnetic diffusivity of metallic hy- 
drogen is ?7 « 10^-10"^ cm^s~^. According to Ohm's law 
and Maxwell's equations (hydromagnetic induction equa- 
tion) , a magnetic field will decay unless a velocity field can 
counteract or balance the diffusive effects. For typical val- 
ues of the convective velocity (see above), the characteristic 
magnetic Reynolds number over a star-size conducting re- 
gion is Rm = Vconv R* I V ^ 100. Accordiug to dynamo 
theory, Rm in LMS and BD interiors is thus large enough 
for dynamo to occur, providing both rotation and convec- 
tion are present. Once the criterion for dynamo onset is 
satisfied, the field will grow and is supposed to equilibrate 
when the Lorentz and the Coriolis forces become compa- 
rable (Elsasser number of order unity), reaching the afore- 
mentioned MAC balance. This yields an amplitude for the 
internal field, B^q « (STrpTyfJ)^/^ > 10 G, for fast rotators 

(vrot > lOkms^^). 

At large i-e. in the dissipationless regime, how- 
ever, magnetic diffusion is mainly due to turbulent rather 
than molecular diffusion, with t] = rjt ^ I UconiEl- In that 
case, non-linear saturation occurs when turbulence, enhanc- 
ing the diffusive processes, is strong enough to reduce i?,„ 
down to the critical value for dynamo action, _/?,„> 50. 
This yields typical magnetic length scales ^<i?^/50, and 
thus B, 



eq 



a few 10* G for LMS average conditions, in 
good agreement with surface field determinations of a few 
kC (jPonati et al. 20061 iReiners fc Basri 2007p . This corre- 
sponds to Alfven velocities va 



so that for 



such strong fields the Lorentz force will impede the convec- 
tion by reducing the flow speed. The main effect of a strong 
magnetic field is to inhibit motions across it in comparison 
with motions along it. An ideally conducting fluid is tied to 
the fluid lines. In a fluid of finite conductivity, motion across 
the field is possible at a rate governed by the conductivity. 
In a highly conductive medium with a strong magnetic field, 
the motion will be slow. Stevenson's (1979) stability anal- 
ysis shows that the combination of fast rotation {Ro<^ 0.1) 
and a magnetic field with finite diffusivity enhances con- 
vection, because of the reduction by the Lorentz force of 
the flow anisotropy due to the Proudman- Taylor constraint. 
Stevenson's approach, however, applies to planar geometry. 



^ Observations of sunspot decay indeed suggest that the solar 
surface diffusivity is many orders of magnitude higher than the 
atomic value. 



i.e. thin convection zones, and to uniform density and mag- 
netic fields and is likely to break down for large (star-size) 
convective zones and strong fields. Magneto-convection 3D 
simulations indeed show that the magnetic field inhibits the 
mag nitude of the veloc ity fluctuations and reduces the heat 
flux (jStein et al. 1992[) . On the basis of the aforementioned 
field strength values, it seems unavoidable to suppose that, 
even though the magnetic field will not necessarily stabilise 
the fluid against convection in a fluid of finite electrical 
resistivity, it will cause a serious reduction of convective 
efficiency. 

These estimates, on the other hand, show that the mag- 
netic pressure in LMS or BD interiors, < 10''dyncm~^, 
is orders of magnitude lower than the gas pressure, 
Pg^ 10^2-1016 dyncm-2 (|CBOO|) . and it can be safely ig- 
nored in the internal structure equations. Indeed, if the fluid 
is convectively unstable, the ratio of the magnetic pressure 
inside the flux tube over the surrounding mean gas pres- 
sure is expected to be of the order of the superadiabatic- 

ity, ^'^^ly < (V - V„rf)<10-7 dCough fc Tayler 1966[ 

Me yer 1994] ). An estimate for complete inhibition of con- 
vection is obtained when the Lorentz force is strong enough 
to balance the buoyancy force: 



AttI 



> pgS(y - Vad) , 



(2) 



where I <C R* is the aforementioned characteristic mag- 
netic length scale in a dynamo turbulent medium. This 
yields field amplitudes of the order of lO'' G, compa- 
rable to the value of Bgq. Magnetic fields can thus in 
principle severely inhibit convection in the interior of ac- 
tive LMS and BDs. The criterion ^ for stability against 
convection in the presence of a magnetic field is simi- 
lar to the one derived by Stevenson in the dissipation- 
less regime and by Gough & Tayler (their Eq.(1.2)), ex- 
cept for the reducing factor ~ //i?* for the gas pressure 
term (assuming VPg ^ Pg/Ri,). The Gough and Tayler 
approach, however, is primarily devoted to the study of 
surface spots, where the magnetic and gas pressures are 
comparable. Applying this criterion to LMS interiors yields 
field strengths > 10^ G for the inhibition of convection in 
the core of a 0.3 Mq star (jMuUan fc MacDonald 2001)) . 
It seems rather difficult to generate such strong fields. 
Together, near-equipartition (within a factor ^10) between 
turbulent and magnetic energy, B^/Stt ^ pVconv^ and the 
fact that L{r) cx Airr'^ (pv^^^^) is a slowly varying quan- 
tity, yield an amplification factor ^10-100 from the sur- 
face to the central regions. A full field of several mega- 
Gauss would thus be in super-equipartition and, if con- 
fined to the interior, would be unstable (Elsasser numbers 
> 1) dTayler 19731 |Markey fc Tayler 1973[ ). Therefore, it 
does not seem realistic to apply the Gough- Tayler criterion 
to the entire stellar structure, in particular for uniformly 
dense objects like LMS and BDs. 

A real picture of a magnetised, convectively unstable 
medium is cooling flows along the magnetic flux tubes pass- 
ing through the convective medium. Part of the thermal 
flux, confined to the tubes, is thus carried by diffusion. This 
contraint on the flow patterns leads to a substantial reduc- 
tion in the transport of energy and reduces the maximum 
possible heat flux. Given the absence of a proper treatment 
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of heat transport in a magnetised medium for LMS con- 
ditions, which requires multidimensional simulations over 
a characteristic convective length scale, and given the ex- 
ploratory nature of this Letter, we decided to stick to a 
minimalist approach, focusing on the reduced efficiency of 
global thermal convection due to the presence of strong ro- 
tation and/or magnetic fields in EBs. Such an approach is 
based on the phenomenological representation of this net- 
work of diffusive tubes surrounded by impeded field-free 
convective regions as a global convective system - since 
the large Reynolds numbers illustrate the overwhelming 
importance of macroscopic motions over microscopic dif- 
fusive processes - with reduced efficiency. In the framework 
of the standard MLT formalism, this translates into a mix- 
ing length parameter a = I /Hp < 1. 



3. Effect of spot coverage 

The strong chromospheric and coronal emission in LMS of 
spectral types > M2-M3, with Lh„ and Lx about 100 times 
the solar value (jGizis et al. 20021 [Mohanty fc Basri 2003] ), 
can be associated with an average large fraction of the radi- 
ating surface being covered with magnetic spots. Cool spots 
are the illustration of the inhibition of energy transport 
by convective motions (buoyancy) in a (rotating) highly 
conducting medium, but they also illustrate the fact that 
the field is not strong enough to suppress convection com- 
pletely over the entire structure. This brings some justifi- 
cation to our admitedly simplistic approach where the aim 
is to explore the effect of reduced convective efficiency in 
EBs. We denote f3 = Ss/S^, as the (time-averaged) fraction 
of the stellar surface covered by cool spots, the total 
flux emerging from the spots, and J-i, = aTes^ the one as- 
sociated with a spot-free star of effective temperature Toff- 
Cool spots imply Ts < (iTi,. The total flux is 



T - (1 - /?) + < T,. 



(3) 



Cool spot coverage thus yields an effective temperature 
Teff < Teff^,, a c onsequence of the reduced heat flow reach- 
ing the surface ( Cowling 1976[IStein et al. 1992p . Given the 
presently undetermined spot temperature for LMS and the 
exploratory nature of the present calculations, we simply 
assumed black (zero-temperature) spots. The effect of spot 
coverage for a given value of (3 is thus an upper limit. 



4. Results 

4.1. Effect of reduced convective efficiency 

Figure [Ha) portrays the m-R relationship for LMS at 1 
Gyi0 for various values of the mixing length parameter a. 
Observed radii of EBs are displayed as well. Objects below 
~ O.35Af0 are fully convective (Chab rier fc Baraffe 1997]) . 
convection is nearly adiabatic, so changing a for these ob- 
jects has a modest impact. Since superadiabaticity increases 
with mass, in particular in the outermost regions of the 
star, reducing the convective efficiency, i.e. decreasing a, 
for higher masses leads to an increasingly greater effect. 
The reduced convective flux requires a larger fraction of 




a = 0.5, /S' = 0.5 
a=\.0, S = 0.3 J 
a = 2.0, ^ = 0.0 f 




Objects with m > O.OQA/© are on the ZAMS at this age 



(fCBOOl) . 



Mass [Mg] 

Fig. 1. Mass-radius relationships for 1 Gyr old LMS for 

various values of the mixing length parameter a=l/Hp ((a) 
upper panel) and of the fractional surface area (3 covered 
with cool spots ((b) lower panel). Inset upper panel: dash- 
dot line : a = 2, 10^ yr. Solid thick line: standard case, 
no spot coverage. Low-mass eclipsing binary observational 
determinations are displayed with their la error bars. 



the heat to be transported by radiation in the outermost 
regions, yielding a steeper outer thermal gradient and thus 
a cooler T^ff. The immediate consequence is to decrease the 
luminosity and thus the central temperatures, i.e. the nu- 
clear energy production needed to maintain thermal equi- 
librium. This in turn yields an expansion of the star, i.e. 
a larger radius. For BDs, the reduced heat flux implies a 
slower contraction rate. 

We have also explored the fate of a 0.3 Mq star when 
drastically reducing the convective flux in the interior (a = 
0.05) as a result of condition 0. The star indeed develops 
a stable inner radiative zone (Vrod < Vad) over ~ 20% of 
its mass, yielding a 7% larger radius and a cooler Toff. 

4.2. Effect of spot coverage 

Figure [IJb) illustrates the effect of surface spot coverage. 
As mentioned above, spot coverage yields a smaller heat 
flux output, i.e. cooler Tgg and thus larger radii (since L is 
unaffected) at a given age compared with spot-free objects. 
A spot coverage fraction (3 ~ 30-50%, for a value of a = 1, 
can by itself reproduce most of the observed EB radii within 
Icr. 

Figure [2] displays the effective temperature as a func- 
tion of the radius for low-mass EBs for different values of a 
and p. Both reduced convective efficiency and spot cover- 
age yield larger radii and cooler Toff compared with regular 
objects and provide a simple explanation for the partic- 
ular mechanical and thermal properties of EBs. Figure [3] 
portrays m-R relations for the recently discovered eclipsing 
BD system (I Stassun et al. 2006^ . A solution with a = 0.5, 
/3 = 0.5 for the most massive object (mi=0.054MQ), 
and a=l, /3=0.2 for the least massive one (m2=O.O34M0) 
yields radius values within the error bars and a temper- 
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Radius [Rg] 

Fig. 2. Effective temperature vs radius at 1 Gyr for various 
values of a and /3. Solid thick line: standard case, no spot 
coverage. 
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Fig. 3. Same as Fig. 1 for the eclipsing BD system (Stassun 
et al. 2006) for different values of a and /3. 

ature ratio Tcffa/Tcff i=2440 K/2320 K=1.05, in agreement 
with the observational determination, for the age of the 
system ~ 10^ yiQ. This possible solution suggests that 
the evolution of the most massive BD has been signifi- 
cantly affected by the presence of a strong magnetic field, 
whereas the least massive BD has been less affected. Indeed, 
BDs with spectral type >M8-M9, i.e. Toff < 2400 K, are 
rapid rotators and are believed to shelter large-scale mag- 
netic fields in spite of showing no sign of persistent activ- 
ity dMohanty fc Basri 2003[ IChabrier fc Kiiker 2006|) . The 
reason invoked for this is the low electrical conductivity 
and thus the lack of substantial current generation in their 
atmosphere ( jMohanty et al. 2002[ IReiners fc Basri 2007|) . 
Spot coverage for cool (Toff < 2400 K) BDs is thus expected 
to be small, even in the presence of a strong magnetic field. 
A possible, although speculative, explanation for the BD 
system of interest is that the most massive object was orig- 
inally hot enough (its Tes for a=2, /?=0 is Teg = 2800 
K) for significant spot coverage to occur, whereas the least 

^ The larger Tcff determinations obtained by Stassun et al. 
(2006) with the usual Sp-Tcft-colour and Sp-gravity relations 
suggest that these relations cannot be applied to active objects 
and to stellar surfaces exhibiting magnetic spots. 



massive BD remained too cool for atmospheric coupling 
between the magnetic field and the gas to be significant, 
yielding less spot coverage. 

5. Conclusion 

In this paper, we have examined the consequences of (i) 
inhibiting convection due to rotation and/or internal mag- 
netic field and (ii) the presence of surface magnetic spot 
coverage, on the evolution of LMS and BDs. We have fo- 
cused on the particular case of EBs. Our approach is phe- 
nomenological and the present paper makes no claim to 
present a consistent description of the effect of a large 
magnetic field on the evolution of dominantly convective 
objects. A proper approach to such a complex process for 
LMS interior conditions requires presently unavailable nu- 
merical tools. Our calculations show that rotation or mag- 
netic field induced inhibition of convection leads to a re- 
duced heat flow and thus significantly larger radii and 
cooler Teff than for regular objects. Spot coverage of the 
stellar radiating surface by itself has a major impact on 
the evolution. Either one or a combination of these effects 
explains the observed mechanical and thermal properties 
of EBs. Work is under way to explore the spot modu- 
lation in EBs and to derive observational values of the 
spot coverage fraction /3 for M-stars (Morales et al., in 
preparation). A spot coverage fraction /3 ^ 30-50%, how- 
ever, is consistent with observations of rapidly rotating, 
more massive active stars (J effers et al. 200"5!|) . These calcu- 
lations strongly suggest that the evolution of EBs, or of very 
magnetically active LMS and BDs, differs noticeably from 
the one of objects for which rotational or magnetic effects 
are negligible, a fact supported by observational analysis 
([Morales et al. 2007P . This should be taken into considera- 
tion when comparing observations with theoretical models. 
Global averaging of the mass-radius relationship for LMS 
is thus incorrect and highly misleading. We show that the 
puzzling eclipsing brown dwarf system that was recently 
discovered, with the more massive companion being the 
cooler one, can be explained if the evolution of the most 
massive object has been strongly affected by magnetic ef- 
fects. 

We suggest that internal field strengths consistent with 
the observed surface value, about a few kG, might be 
sufficient for severely inhibiting convection in parts of 
LMS interiors, pushing the limit for the onset of an in- 
ner radiative zone below ~ 0.35Mq, with important con- 
sequences for the field geometry ([Chabrier fc Kiiker 20061 
IDonati et al. 2006p . 
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